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We study the problem of spatially stabilising four dimensional extremal black holes in background 
electric/magnetic fields. Whilst looking for stationary stable solutions describing black holes kept in 
external fields we find that taking a continuum limit of Denef et al's multi-center solutions provides 
a supergravity description of such backgrounds within which a black hole can be trapped in a given 
volume. This is realised by levitating a black hole over a magnetic dipole base. We comment on 
how such a construction resembles a mechanical Levitron™. 



I. INTRODUCTION 

Being motivated by on-going interest in questions con- 
cerning black hole production; in this note we address 
a curiosity regarding how one could go about stabilising 
such a black hole using external fields, thus leading to a 
black hole analog of a particle-trapor rather as we shall 
see that of a Levitron™. However unlike the more fa- 
miliar subatomic particle traps or even Millikan's famous 
oil drop experiment [T], the effects of general relativity 
give rise to interesting new features. We shall describe 
how this idea can in fact be materialised by writing down 
solutions for black holes levitating in electromagnetic as 
well as constant gravitational fields. 

For the purpose of this note we consider four di- 
mensional extremal black holes with electric/magnetic 
charges q and p respectively. Extremal black holes are 
BPS solutions to four dimensional supergravity. The 
most general metric ansatz consistent with supersymme- 
try can be written as ( [2], [3], [5] ) 
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is the four dimensional gauge field. 'P(x), Q(x) are har- 
monic functions associated to charges p and q respec- 
tively. is the Dirac part of the vector potential satis- 
fying d<d = *dV(x) with the Hodge star * defined on ]R 3 . 
For a single spherically symmetric black hole in vacuum, 
it holds that 3 = 0. However for our considerations, 
we shall be looking for solutions when the black hole is 
placed in external electric and magnetic fields. There is 
now a non-zero Poynting vector corresponding to a ro- 
tating geometry. We first look for levitating solutions in 
constant background fields. It turns out these are inad- 
equate for stabilisation in all three directions. Then we 
look for more non-trivial backgrounds obtained using a 
continuum limit of Denef et al's [2J, [3J, [5] multi-center 
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solutions and find that turning on dipole fields achieves 
the desired result. 



II. BLACK HOLE LEVITATION IN CONSTANT 
EXTERNAL FIELDS 

Given the metric ansatz in eq.fJTJ, we begin by looking 
for stationary solutions of a black hole placed in con- 
stant electric, magnetic and gravitational fields. In order 
to achieve this we have to specify explicit harmonic func- 
tions describing this configuration, then compute the off- 
diagonal elements lj and solve the associated integrability 
equations. We claim that the desired harmonic functions 
describing this configuration are 



V(x) = 
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where B and E are constant magnetic respectively elec- 
tric fields oriented along the z-direction and z denotes 
the z-coordinate. I marks the position of the black hole's 
horizon, which we determine via integrability conditions. 
u, v are constants. In principle, we can absorb u and v 
via a shift in the z-coordinate. This point will be made 
clear when we solve for I. The Bz and Ez in eq.fj^ 
are linear terms that satisfy Laplace's equation and can 
be recognised as the usual electro/magneto-static poten- 
tials associated to constant fields. Note that extremality 
implies the above linear terms also source constant grav- 
itational fields. 

A nice way to motivate the expressions for V(x) and 
Q(x) is to extract them via a special limit of Denef et 
al's multi-center solutions [2], [3J, [5]. More specifically, 
let us consider the two-center solution. This is a regular 
BPS solution of four dimensional Af = 2 supergravity. It 
is stationary but non-static and hence caries an intrinsic 
angular momentum. Moreover the black holes compris- 
ing this bound state possess mutually non-local charges. 
Let us denote the corresponding two charge vectors as 
r = (p, q) and r = (p, q). The idea is now to carry 
the charge T all the way to infinity while scaling (p, q) 
and the radial coordinate of the charges in such a way 
that the magnitudes of the electric/magnetic fields them- 
selves are held fixed. Applying this limit to the expres- 
sions for electro/magneto-static fields of point charges 
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indeed leaves us with constant fields oriented opposite to 
the direction of the source charges T. Without loss of 
generality, the z-axis can then be chosen to point in the 
direction of the sources. Integrating these fields along the 
line element, precisely yields the linear potential terms in 
eq.@. 

In fact we may also use this limiting two-center system 
to captures other features of our original configuration of 
a black hole in constant external fields. Following 0, [5], 
[5], we can determine the off-diagonal terms in the metric 
using 



V x uj = V{x)VQ(x) - Q{x)VV(x) 



(3) 



Below we shall solve to for a class of precessing solutions. 
Furthermore operating a gradient on both sides of eq. ^ 
leads to the following integrability equation 



V(x)V 2 Q(x) - Q(x)V 2 V(x) = 
which we evaluate at x = I to get 



qu — pv 
pE-qB 



(4) 



(5) 



This gives us the position of the black hole. Here I = 
(0, 0, 1) can be chosen on grounds of symmetry. One can 
also perform a shift of coordinates, so as to place the 
black hole at the origin. This can be achieved by setting 
constants u = v = 0. Note however that (pE — qB) ^ 
is required in order to preserve mutual nonlocality. 

Eq.([3]) can be conveniently solved using spherical co- 
ordinates (r, 9, <f>). And that leads to a system of coupled 
differential equations 
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while (V x Q) , = due to ^-independence on the right- 
hand side. Our objective is now to seek out a non- 
trivial solution which confers to the description of a black 
hole rotating in the presence of external electromagnetic 
fields. We find that there exists such a simple solution 
with azimuthal symmetry 



= sin 9 (pE — qB) 



(T) 



while uj r = tog = 0. For completeness let us also mention 
that the solution presented in eq.([7]) is certainly not the 
most general. For instance, we also find that solutions 
with harmonic variations such as = cos <p also exist 
and very likely one may well find a more general class of 
these. But we shall not require that for our purposes. 

The solution above allows us to levitate a black hole at 
a fixed height on the sy-plane owing to the balancing act 
between gravitational attraction and electro/magneto- 
static repulsion. However it is not stable in all three 
directions and can move about the surface of the plane. 
To localise the black hole in all three directions we need a 
more complicated background field where the black hole 
can be held at a local minimum of an effective potential. 



III. CONTINUUM LIMIT OF MULTI-CENTER 
SOLUTIONS 

In this section we start looking for extremal stationary 
solutions to Einstein-Maxwell gravity that admit back- 
grounds with multipolc electromagnetic fields. As be- 
fore, we work with four dimensional gravity with just 
one gauge field. Generalisations to n — 1 vector fields or 
inclusion of other charges such as DO and/or D6 in Type 
II A are rather straightforward. Let us now see how 
taking a continuum limit of Denef et al's multi-center so- 
lutions yields the desired backgrounds. In order to write 
down harmonic functions for such a smeared distribution 
of black holes, we define density functions p e [x'), p m {x') 
via 



p e (x')d,T' 



and 



p m (x')dr' = P (8) 



where dr' is a volume element within a compact support 
V, that covers the distribution. In the continuum limit, 
harmonic functions for multiple black holes take the form 



Q{x) 



Pe{x') 



dr' P{x*) = u + 



Pm{x') 



dr' 
(9) 

To these harmonics one may also add linear terms Ez and 
Bz corresponding to constant fields, whenever required. 
From a computational point of view, the real utility of 
the above-mentioned smeared distributions shows up in 
their respective multipole expansions. Expressing this in 
the regime that \x\ » \x'\ holds, we have 



Q(x) = v 
V{x) = u 
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where Q, P are electric respectively magnetic monopole 
moments; A e , A m are electric and magnetic dipole mo- 
ment vectors; and T e , T m are respectively electric and 
magnetic quadrupolc moment tensors - all defined in the 
usual way. We employ boldface characters to denote vec- 
tors as well as tensors. The " " in eq.(10l denote 



terms with higher order moments. When \x\ » \x'\, the 
series is convergent and these functions can be used to 
describe supergravity solutions associated to any specific 
multi-moment source, provided all lower moments van- 
ish for that distribution. As an illustrative example, we 
analyse the solution for a charge distribution with dipole 
order corrections. 

First let us check that the functions in eq. ^ yield 
integrability conditions that are well-defined in the con- 
tinuum limit. Evaluating eq.Q for these harmonics gives 



Pe (x)T(x)- Pm (x)Q(x) = 



(11) 



Outside the support V, this expression vanishes identi- 
cally; whereas points x € V inside the support region 
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yield 



IV. TOWARDS A BLACK HOLE LEVITRON™ 



Up e (x) - Vp m (x) = 



(12) 



In order to compare this to the analogous result for a 
single black hole, we integrate both sides of eq.(12l over 
the volume V to get 



uQ - vP = 



(13) 



This is exactly what one has for a single-center solution 
with charges Q and P; thereby confirming the asymptotic 
dependence of u and v for an arbitrary multi-center con- 
figuration having fixed total (monopole) charges Q and 
P. 

Having checked consistency of integrability conditions, 
we next compute the off-diagonal elements Co in the met- 
ric via 



V x a; 



-V(x)E(x) + Q(x)B(x) 



(14) 



where E(x) and B(af) refer to exact electric and magnetic 
fields corresponding to distributions p e (x) and p m (x) re- 
spectively. In this sense the continuum limit described 
here is much simpler than a finite N many body black 
hole system for which integrability equations turn out to 
be quite hard to solve in full generality. 



For our objectives, it will suffice to solve eq.(14l using 



its multipolc expansion. As an illustration, we consider a 
smeared distribution where the monopole contributions 
to Co get magnetic dipole corrections coming from A m , 
which is aligned along the z-axis. In spherical coordi- 
nates, eq.(14l takes the form 



(V x Q) r 
(V x uj) e 



2 v A m cos 9 ^ Q A m cos 9 
v A m s'm9 Q A m sin# 



(15) 



while (V x lu)^ — due to symmetry in the </>-direction. 
Note that whilst writing down eq.(fT5| , we make use 
of the integrability constraint eq.(|13[) (inserting it into 
eq.(14| ). As before, we seek solutions characterised by 



azimuthal symmetry. The ensuing result is 
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(16) 



and ui r = tog = 0. At large distances away from the 
smeared sources, eq.(16) gives dipole corrections to lead- 



ing order contributions in the metric. In fact these con- 
stitute sub-leading contributions to the geometry. It is 
these multipole corrections that distinguish a true one- 
centered black hole from a multi-center distribution of 
black holes, when viewed at asymptotic infinity. For a 
pure one-center solution, Co identically vanishes. While 
for the multi-center case, it is non-trivial but quite dif- 
ficult to compute for any given discrete configuration. 
The continuum limit, on the other hand, facilitates vi- 
able computations, at least order by order in a multipolc 
series expansion. 



We are now ready to combine results of the last two 
sections to construct stable levitating black hole solu- 
tions and realise a Levitron™-like construction. We per- 
turb the constant background fields of section II with a 
magnetic dipole field and over this perturbed background 
solve for a black hole held at a fixed height. The dipole 
fields are produced by the smeared distribution discussed 
in section III. For simplicity we consider a black hole with 
only electric charge q ( a dyonic generalisation is also 
straightforward ). This construction is captured by the 
following harmonics 



Q(x) = v 
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Ez V(x) 



Am COS ( 
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The dipole moment is aligned parallel to the z-axis and 
carries a magnitude A m . While 9 is a coordinate denot- 
ing the angle that the position vector x makes with the 
z-axis. However it will suffice to turn off the constant 
fields E and B for the rest of the computation. As we 
shall see this is because a dipole background is sufficient 
to hold the black hole at a fixed height and keep it sta- 
ble in all three directions. Superposing constant fields 
do not affect stability of the solution but ultimately we 
will need the constant fields for giving an interpretation 
of black hole levitation in a constant gravitational field ( 
as would be the case if we were ever to trap a small black 
hole in a laboratory somewhere on Earth ! ). 

Continuing with the calculation, the position of the 
black hole I is determined by evaluating eq.jijl at the 
location of the pole x = I using harmonics in eq.(17l 
with E = B = 0. This gives 



M = 



— A m cos( 
it 



(18) 



Also evaluating the integrability equation at the other 
pole x = determines the constant v 



(19) 



Physical solutions only exist for I (= | Z |) real and non- 
negative. For instance, when u > 0, then 9 can at- 
tain values from to -| provided the dipole is directed 
along the negative z-axis. <j) remains unconstrained. On 
the other hand, for a dipole pointing in the positive z- 
direction, the angle 9 spans within the range ^ to tt. 
(see fig. 1 below). When u < the signs appropri- 
ately reverse. The solution space of the black hole is now 
confined to a restricted parameter space. More precisely 
these are circular orbits corresponding to given values of 
9 on an equipotential surface of a dipole field. And in 
turn each orbit refers to a solution with a specified radial 
distance I. 

In fig. 1 below, we plot eq.(18l. At 9 = the black 



hole sits at a fixed height on the z-axis; at 9 = \ it falls 
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into the origin; while the case < < | corresponds 
to the black hole being located anywhere on a circular 
orbit centered at height I cos and having radius I sin 9. 
Solutions on the positive z-axis correspond to the case 
when A m < (for u > 0), while those on the negative 
axis refer to A m > 0. For each value of 9 in eq.( 18 1 there 
exists a solution for uj. At 9 — the solution space is just 
a single point and that is when the black hole achieves 
stability in all three directions at a fixed height on the 
z-axis. 

For completeness we first compute Q when the black 



hole is still sitting at the origin, that is when / = 0. After 
that we shall determine the modification in Co required to 
achieve stable levitation at a fixed height on the z-axis. In 
fact the solution at I — 0. can simply be borrowed from 
our calculation in eq.( 16 ) once we make the substitutions 
Q -> q and P -> 0. 



On the other hand, when the black hole is made to 
levitate at a fixed height I on the z-axis we have to solve 
the following system of equations 
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(Vxw) e = 
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(20) 



and again (V x uj)^ = 0. Also I = (0,0,1). This now 
becomes fairly more complicated compared to the non- 
levitating case. The modification in the metric reflects a 
modification to the geometry of the system. If we restrict 
to azimuthally symmetric cases, we find that eq.(20l has 
a solution only for small heights of levitation, that is 
when I « r. This can be understood in the following 
way. In this set-up the system consists of the black hole 
plus the source of the dipole field. Let us call the latter 
the base. The levitating we arc looking for requires that 
the base be rigid against the gravitational pull of the 



black hole, that is the center of mass of the whole system 
be as close to the base as possible. For very large charges, 
corresponding to large values of I, a stable symmetric 
levitating solution does not seem to exist ( we see this 
from numerical checks ). In that case more complicated 
non-symmetric solutions may be sought for, but we would 
hardly call those levitating. 

Narrowing down to our regime of interest, we expand 
around I « r and solve eq.(20l order by order in I. 



Truncating up to second order terms we get 
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while uj r = tug = 0. This solution enables us to write 
down the full metric for a stationary system of a black 
hole levitating in equilibrium above a magnetic dipole 
field. Also this calculation easily extends to the case of 
a dyonic black hole. 

Comparision to a Levitron™ 

We now compare the levitation of black holes discussed 
above with that of a Levitron™ The latter is a 
spin stabilised magnetic levitation device first invented 
by Roy Harrigan 6 . It basically consists of a permanent 



base magnet above which a spinning top with a magnetic 
dipole moment levitates mid-air and is stable in all three 
directions. This gives rise to an apparent paradox due to 
Earnshaw's theorem [H] which states that no stationary 
configuration composed of electric/magnetic charges and 
masses can be held in stable equilibrium purely by static 
forces. And the reason for this is simply that all static 
potentials satisfy Laplace's equation whose solutions only 
exhibit saddles at critical points : there are neither any 
maxima nor minima. It was Sir Michael Berry's [9] (see 
also [TU]) remarkable insight invoking adiabatic averag- 
ing that helped resolve the apparent paradox. He showed 
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that a slow precession mode ( when averaged over the fast 




FIG. 1: Plotting the solution space for / 

rotation mode ) was responsible for creating an effective 
stationary potential with a stable minimum. This is the 
same principle used in neutron traps as well as other par- 
ticles carrying magnetic dipole moment. 

A natural question which arises is whether our 
black hole construction also mimics the physics of the 
Levitron™ and how it overcomes Earnshaw's theo- 
rem. The latter it already seems to evade since it is 
based on Einstein's gravity rather than Newton's. How- 
ever the gravitational interpretation of our Black Hole 
Levitron' s™ balancing mechanism admittedly requires 
further investigation. Nevertheless a naive classical intu- 
ition can be obtained from the fact that a non-vanishing 
Poynting vector gives rise to a rotating black hole geom- 
etry and in turn a rotating electric distribution induces a 



magnetic field that repels the base magnet. It is the uj in 
the metric that is responsible for inducing this balancing 
force. On the other hand the gauge theoretic interpre- 
tation of this multi-center balancing has been better un- 
derstood in terms of Denef 's quiver quantum mechanics 
[4] wherein the distance between centers is determined 
via an effective potential whose minima determine the 
stability loci I. 



V. CONCLUSIONS AND DISCUSSION 

In this article we have constructed a levitating black 
hole solution. Our Black Hole Levitron™ stabilises an 
extremal black hole at a fixed location in an electromag- 
netic field produced by a continuous distribution. Our 
work is built-up using Denef et al's multi-center solutions, 
which by themselves arc stable, stationary BPS solutions 
with non-local charges. Our harmonic functions and in- 
tegrability conditions can all be retrieved as special limits 
of the discrete multi-center case. Therefore our levitat- 
ing solutions also describe stable, stationary configura- 
tions. This black hole construction very much resembles 
a mechanical Levitron™ and it would be interesting to 
investigate if Berry's mechanism can be proven to apply 
to this set-up as well. And finally it would be of prac- 
tical relevance ( in future ! ) to construct solutions for 
non-extremal Black Hole Levitrons™ ! 

Other interesting directions might be further investi- 
gation into other applications of the continuum limit of 
multi-center solutions. Compared to discrete-centered 
configurations, the smeared distribution lends itself to 
more viable computations. One may ask what role these 
distributions play in microstate counting of multiple- 
black hole geometries. 
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